Helical magnetic state in the distorted triangular lattice of 

a-CaCr204 



L.C. Chapon/ P. Manuel/ F. Damay,^ P. Toledano,^ V. Hardy,^ and C. Martin^ 

^ISIS facility, STFC Rutherford Appleton Laboratory, 
Chilton, Didcot, Oxfordshire, 0X1 1-OQX, United Kingdom 
Lahoratoire Leon Brillouin, CEA-CNRS, UMR 12, 
CEA-Saclay, 91191 Gif-sur-Yvette Cedex, France 
^Laboratory of Physics of Complex Systems, 
University of Picardie, 33 rue Saint-Leu, 80000 Amiens, France 
^Laboratoire CRISMAT, ENSICAEN, UMR 6508 CNRS, 
6 Boulevard du Marechal Juin, 14050 Caen Cedex, France 
(Dated: June 30, 2010) 

Abstract 

The magnetic properties of tlie high temperature a form of the CaCr204 compound have been 
investigated for the first time by magnetic susceptibihty, specific heat and powder neutron diffrac- 
tion. The system undergoes a unique magnetic phase transition at 43K to a long range order 
incommensurate helical phase with magnetic propagation vector k = (0, 0.3317(2), 0). The mag- 
netic model proposed from neutron diffraction data shows that the plane of rotation of the spins is 
perpendicular to the wave- vector, and that the magnetic modulation is consistent with two modes 
belonging to distinct irreducible representations of the group. The magnetic point group 2221' is 
not compatible with ferroelectricity unlike the CuCr02 delafossite [Kimura et al, Phys. Rev. B, 
78 140401 (2008)] but predicts the existence of quadratic magnetoelectric effects, discussed based 
on a Landau analysis. 

PACS numbers: 71.30.+h, 75.47.Lx, 61.12.-q, 75.10.-b 
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I. INTRODUCTION 



In recent years, spin-driven ferroelectricity given rise to giant magnetoelectric effects 
has been found in a variety of frustrated magnets, in particular in systems that develop 
long-wavelength magnetic modulations such as cycloidal or helicoidal states, incommen- 
surate or not with the crystalline lattice. Most of the attention was initially focussed on 
cycloidal magnets such as TbMnOs-^ for which the inverse Dzyaloshinskii-Moriya^ and spin- 
current modeF predict the onset of a spontaneous polarization when the direction of the 
modulation (k-vector) is perpendicular to the spin-rotation axis (e)^. More recent work on 
incommensurate magnetic phases in triangular lattices^"— has established that improper 
ferroelectric order of similar magnitude also appears in spirals magnets, i.e systems that 
possess a true magnetic chirality (k parallel to e), however requiring another microscopic 
coupling mechanism such as hybridization effects proposed by ArimaM. Those systems are 
quasi two-dimensional triangular lattices and display a nearly 120° magnetic arrangement 
in the triangular layers as expected for Heisenberg spins. Improper ferroelectricity predicted 
by symmetryl° arguments depends only on the orientation of e with respect to the crystal 
axes, often dictated by weak single-ion anisotropy terms. 

In this paper, we report for the first time on the magnetic properties of a;-CaCr204 
whose orthorhombic crystal structure^^ is closely related to the structure of delafossites, 
at the exception of a small distortion of the triangular lattice. The magnetic suscepti- 
bility indicates strong antiferromagnetic correlations (Weiss temperature of -893K) and a 
unique transition to a long-range ordered state at TAr=43K. The specific heat displays a 
sharp and intense peak around 43K, consistent with a unique transition. Powder neutron 
diffraction proves the existence of long-range magnetic order below Tn, incommensurate 
with the crystal lattice with propagation vector k=(0,0. 3317(2), 0). The model proposed 
for the magnetic structure, derived from the neutron data and symmetry considerations, 
corresponds to a nearly 120 ° out of plane helix, with e perpendicular to k. Symmetry 
analysis using the complete irreducible co-representations of the wave-vector group shows 
that the point group symmetry of the magnetic state is 2221' and can be only stabilized 
through two successive second-order transitions or through a first-order transition. Whilst 
in the rhombohedral delafossite such out-of-plane spiral leaves a polar point group, here the 
spin chirality preserves all proper rotations, predicting the absence of electric polarization 
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at Tat. However, we show through a Landau analysis of the magnetoelectric free energy 
that quadratic magnetoelectric effects are allowed by symmetry and should be observed in 
low magnetic fields. 



II. EXPERIMENTAL 

The synthesis of polycrystalline a-CaCr204 followed a two-steps process. A sample of 
/3-CaCr204 was first prepared in the shape of a rod (6 millimeters in diameter and several 
centimeters in length), starting from a 1:1 mixture of CaO and Cr203 heated at 1400°C for 
12 hours in an argon flow. This rod was then heated in argon atmosphere by using an image 
furnace to obtain the high temperature conditions (T>1700°C) but without reaching the 
melting temperature. Magnetic susceptibility has been measured using a Superconducting 
Quantum Interference Device (SQUID, Quantum Design @) in a magnetic field of 0.3T on 
warming following a zero- field cooling process. Heat capacity measurements (C(T)) were 
carried out in a Physical Property Measurement System (PPMS, Quantum Design @) using 
a semi-adiabatic relaxation method. Two types of procedures were followed to record zero- 
field C(T) curves: (i) equally spaced data points registered upon warming and using a small 
temperature rise of the order of 0.4 K around the transition. The analysis is the standard 
one in which the so-called 2r model^^ is used to fit at once both the heating and cooling 
branches at each point; (ii) a large single-pulse method described in^^ where the temperature 
is swept across the complete width of the transition, either warming or cooling (temperature 
rise of 5K). It this case, the heat capacity is derived from a point-by-point analysis of the 
time response along each of the two branches. In case of a first-order transition, it is known 
that the former method can yield spurious results,— whereas the latter is able to account 
for both the latent heat and hysteretical features.— Neutron diffraction experiments were 
collected on the WISH diffractometer of the ISIS Facility, Rutherford Appleton Laboratory 
(UK). Data have been focus on five histograms from detector banks covering 32 ° each, 
and are shown for the bank centred at 2^=90°. A 0.55g sample was placed in a thin wall 
(30/im) vanadium cylindrical can mounted in an helium cryostat. Diffractograms have been 
recorded in the paramagnetic regime at 70K, at 1.5K and warming in varying temperature 
steps. Rietveld refinements of the neutron data have been performed with the FuUProf 
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programis. 



III. RESULTS 




FIG. 1. (Color online) Left) Crystal structure of a-CaCr204. CrOg octahedra are shown in blue 
and green colors for the inequivalent Cri and Cr2 sites (see text for details). The Ca^"*" ions 
are shown as light grey spheres. The red thin line marks the crystallographic unit-cell. Right) 
Projection of the Cr02 layer onto the 6c-plane. Colored thick lines shows the connectivity between 
first-neighbor Cr ions. The corresponding interactomic distances are shown. 

The crystal structure of a-CaCr204 was first reported by Pausch et al.— . It crystal- 
lizes with orthorhombic symmetry, confirmed by our Electron Microscopy measurements 
(not shown) in the space group Pmmn (Fig. [1]). The structure refined from our neutron 
diffraction data in the paramagnetic regime at 70K is also consistent with the earlier work—. 
There are two inequivalent Cr sites, Cri located at (0,0,0), of 1 symmetry and Cr2 at position 
(0. 5047(6), i,0.497(l)) located on a mirror plane, each site of multiplicity four. The Cr ions 
are six-coordinated by oxygen in a distorted octahedral configuration. The octahedra share 
edges to form Cr02 layers in the 6c-plane, leaving a weakly distorted triangular lattice as 
shown by the different interatomic distances marked in Fig. [TJ The Cr02 layers, separated 
by Ca^"^ ions, are stacked along the a-axis of the orthorhombic unit-cell. The in-plane mag- 
netic interactions between Cr can be mediated either by direct exchange through the t2g 
orbitals or by Cr-O-Cr super-exchange interactions. The interactions between layers can be 
mediated by super-superexchange interactions through two oxygen atoms, but are expected 
to be much weaker due to the large interatomic Fe-Fe distance (5.52 A at 1.5K). 

The inverse magnetic susceptibility (Fig. |2]) shows a quasi-linear regime above 250K. A 
linear fit of the data with a Curie- Weiss law in the temperature range 250K to 350K yields 
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FIG. 2. Inverse molar magnetic susceptibility of a-CaCr204 as a function of temperature under a 
magnetic field of 0.3T. The line is a linear fit with a Curie- Weiss law. 

a paramagnetic moment of 4.3/iB per Cr^"^ ion, and a Weiss temperature of -893K. The 
value of the paramagnetic moment is larger than the expected value of 3.87 fiB for a pure 
spin contribution (S=3/2) of Cr'^'*'. The reason for this discrepancy is likely related to the 
limited temperature range of the measurement, much lower than the Weiss temperature, 
and it seems rather unphysical to ascribe this deviation to a large orbital contribution. The 
strongly negative Weiss temperature indicates large antiferromagnetic (AFM) interactions. 
On cooling below 200K,one observe a deviation of the inverse susceptibility from the nearly 
linear regime followed at Tiv=43K by a more abrupt drop indicative of long-range AFM 
ordering. There are no evidence of additional transitions from susceptibility measurements 
at different fields (not shown). The large difference between the value of the Weiss tempera- 
ture and Tjv is inherent to the reduction of dimensionality in this pseudo-2D system as well 
as magnetic frustration imposed by the triangular geometry. 

The left panel of Fig. [3] displays the specific heat, C(T), recorded from the standard and 
single-pulse methods (see experimental section), while the inset is an enlargement exhibiting 
the shift between the warming and cooling branches derived from the latter technique. At 
TAr=43K, a very sharp peak is observed with a full width at half maximum of about 0.5 
K. Using a Debye function to fit the high temperature data as an approximation of the 
lattice specific heat, one can derive the temperature dependence of the magnetic entropy 
(Sm) shown in right panel of Fig. [31 There is a clear kink at Tjv on the Sm(T) curve, but 
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FIG. 3. Left: Specific heat of a-CaCr204 in zero magnetic field measured by the standard (filled 
diamonds), single pulse method (open circles) and SPM cooling (open squares) (see text for details 
of the various techniques). Right: Magnetic entropy obtained by subtracting the lattice contribu- 
tion estimated by a Debye function to the total specific heat. The inset focus on the temperature 
region close to T^v. The dashed line indicates the theoretical limit of 2Rln(4). 

judging the nature of the transition from this feature only is difficult. The Sm(T) around the 
transition (Fig. [3] inset) can indeed be regarded either as a broadened jump expected for a 
first-order transition (FOT) in real materials, i.e. non ideals with the presence of defects- or 
as a knee typical of a second-order transition (SOT). Two other experimental observations 
must be considered: first, the nearly perfect superimposition of the C(T) curves derived 
from the single-pulse and standard method, whilst a substantial part of the latent heat -if 
present- is supposed to be invisible in the latter technique.— Secondly, the temperature shift 
observed between the heating and cooling data in the single-pulse method is very small, and 
can not even be safely ascribed to a genuine hysteresis as encountered in most of the FOTs. 
In fact, it appears that this hysteresis quantitatively corresponds to the expected temper- 
ature lag between the sample and the platform which contains the thermometer. At this 
stage, one must conclude that the transition at T^r is either a SOT or a very weekly FOT. 
Further information might be derived from the analysis of the critical behavior around T^r, 
but this will require measuring C(T) for an isostructural non- magnetic compound to provide 
an accurate estimate of the lattice contribution more precise than a simple extrapolation of 
a Debye function. This work is currently in progress. 

Below 43K, long-range magnetic order is detected by the presence of five additional 
Bragg peaks in the powder neutron diffraction data, the most intense close to a d-spacing of 
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4.5 A. The intensity of these peaks decreases with the modulus of the scattering vector (Q) 
in agreement with the expected Q-dependence of the magnetic form factor. To determine the 
periodicity of the magnetic structure, an automatic indexing procedure using a grid search^^ 
was employed but failed to give a satisfactory solution due to the presence of pseudo-trigonal 
symmetry. However, the magnetic peaks could all be indexed by exploring reciprocal space 
in the vicinity of the special k-vector expected for a 120° structure in a isotropic triangular 
lattice, corresponding to k=(0,|,0) in the a;-CaCr204 orthorhombic unit-cell. Further analy- 
sis shows that the magnetic structure is actually incommensurate with a propagation vector 
k=(0,0. 3317(2), 0), labelled kg in Kovalev's notation^^. The incommensuration, albeit small, 
is genuine as indicated by the inferior refinement presented in Fig. H] showing an offset for 
the position of the magnetic peak (most noticeable in the difference curve) when the value 
is locked at k=(0,|,0). The profile of the magnetic peaks are almost purely Lorentzian and 
their widths are much larger than the instrumental resolution, indicating short correlations 
lengths. In order to account for the observed profile and derive an accurate value for the 
moment magnitudes, the Rietveld refinement of the magnetic phase was conducted using a 
phenomenological description of the peak broadening by a spherical harmonics expansion 
consistent with the Laue class mmm. 
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TABLE I. Matrix representatives of the complete irreducible corepresentations Ai and A2 asso- 
ciated with the irreducible representations of the Uttle group ri and r2^2 (gee text for details). 
The matrices act each on the two-dimensional spaces spans by the complex order parameters, 
respectively (?7i,??i) and (r/2,^2)- The matrix for time-reversal is not shown. 



Symmetry analysis using the group of the propagation vector (little group) show that the 
magnetic representation for each Cr sites is decomposed in four irreducible representations, 
labelled ri to according to Kovalev's tables^. The magnetic structure is only compatible 
with a model that mix a basis vector along the (1,0,0) direction that transforms as Ti , and a 
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FIG. 4. Rietveld refinements of the neutron diffraction data collected in the paramagnetic phase 
at 70K (a) and in the magnetically ordered phase at 1.5K (b). The data are shown as grey points 
and the result of the refinements as continuous black lines. The curve shifted below represent the 
differences between observed and calculated patterns. In both panels, the top two rows of tick 
marks indicate the positions of Bragg peaks for the phases a-CaCr204 and the impurity Cr203. 
In panel b, the third row of tick marks indicated the positions of magnetic Bragg peaks with 
propagation vector k=(0,0.33,0). 



basis vector along (0,0,1) that transform as T2. In this orthorhombic symmetry, the intensity 
of the diffraction pattern is not sensitive to the relative phase between these two modes, 
i.e. one can not directly discriminate between a transversally polarized spin-density wave 
with the moments pointing in the ac-plane (mixing the two basis vectors with coefficients of 
same characters) or a helical modulation with the moment rotating in the ac-plane (mixing 
the two basis vectors with coefficients of different characters, real and imaginary). However 
at 1.5K, the refined amplitudes of the magnetic modes along the a and c axis were found 
to be 2.31(1) hb- If these modes were in phase (spin-density wave) the amplitude of the 
modulation will largely exceed the fully ordered moment oi 3 /ib expected for a pure spin- 
state Cr3+ ion. On this ground, it appears that the helical modulation is the only possible 
model to account for the observed pattern. It is also physically sensible as the entropy of 
an helical structure with constant moments is more favorable at low temperature than an 
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amplitude modulated structure. The magnitudes of the magnetic moments were found to 
be similar for the two symmetry inequivalent Cr sites, and final refinements were conducted 
imposing such constraint, even if the latter is not directly imposed by symmetry. Finally, 
the relative phase between the modulation for site 1 and site 2 is a free parameter, as these 
sites are not related by symmetry operations of the group. The phase has been refined to a 
value very close to -| and was constrained to this specific value in the final refinement. The 
magnetic structure, showed in Fig. [T|, is made of nearly 120° configuration in each triangle, 
propagating as spirals along the 6-axis of the unit-cell with the two spirals per unit-cell 
(on different Cr sites) of same chirality. Unlike cycloidal modulations, this proper spiral 
structure possess a true handedness. The coupling between adjacent Cr02 layers, mediated 
through super-superexchange interactions, is antiferromagnetic. The magnetic structure is 
similar to that found in delafossites^i^^, i.e. the 120° configuration is out-of-plane. How- 
ever, symmetry considerations lead to different coupling between magnetic order parameters 
and electric order than the delafossite, that will be discussed later. The ordered magnetic 
moment varies smoothly with temperature (Fig. [5]) but the critical exponent was found 
to be ~ ^, deviating largely from the mean- field limit and reminiscent to that found in 
first-order phase transitions^''. Within the temperature resolution of our experiment there is 
a unique transition at Tjv, in agreement with measurements of the magnetic susceptibility 
and specific heat. 

The magnetic model consists of a mixing of two modes belonging to two different irreducible 
representations, i.e. the transition from the paramagnetic phase to the incommensurate 
magnetic phase can be fully described by considering two two-dimensional complex order 
parameters rii=pie'^'^^ ,ril=pie~^'^^ and t]2= P2G^'^^ ,f]2— P2^~^^^ ■ The matrix representatives of 
the complete irreducible co-representations of the group of k are presented in Table HI The 
time-reversal operator 1', which matrix is not shown, needs also to be considered to be 
able to write the Landau energy. Considering these two magnetic order parameters, the 
expansion of the Landau free energy F is written: 



+lplplcos{2(t)) + ... 



(1) 



(2) 
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FIG. 5. (Color online) Left: Magnetic structure of a-CaCr204The two independent Cr sites are 
shown with different colors (light green for Cri and light blue for Cr2). Two unit-cells along the 
6-axis, marked by thin black lines, are shown. The axis are labelled in italic and the direction of 
the magnetic propagation vector (k) is shown by the arrow. Right: Temperature-dependence of 
the ordered magnetic moment on the Cr ions 

, where the a,, /3j (1=1,2) and 7 are the usual coefficients and = 0i — 02 is the phase 
difference between the two modes. Minimization of F with respect to pi and p2 leads to 
five distinct magnetic states that can be stabilized from the paramagnetic space group 
PmmnV . The point groups of these five magnetic states are shown in Fig. [6] together 
with their stability conditions. The structure observed experimentally corresponds to both 
order parameters being non-zero and = 7r/2 (mixing of two modes in phase quadrature) 
and has the point symmetry 2221'. The 1' operator is contained in the point group of 
any incommensurate structure, since its application (phase of vr) is always equivalent to a 
translation. The phase 2221' can be stabilized from the paramagnetic state either through 
two successive second-order phase transitions as indicated by the arrow in Fig. [6] or through 
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a single first-order phase transition. According to our data showing a unique critical point 
but no pronounced hysteresis, we assumed that the transition must be weakly first-order. 
One should note that the only possibility to stabilize the observed magnetic phase through a 
second-order transition would be to consider a symmetry for the exchange Hamiltonian that 
is higher than the crystal symmetry (for example considering an Heisenberg Hamiltonian). 
In such case, the symmetry analysis shows that the two irreducible representations involved, 
Ti and T2, belongs to the same exchange multiplet. However, this is rarely observed and 
always a case of pseudo-symmetry. In the following discussion, we will assume that the 
crystal symmetry is the relevant basis of work, and that the transition is weakly first-order. 

The Landau theory of transitions with two order parameters that are coupled bi- 
quadratically has been studied by several authors^ii^ and reviewed by Toledano et al.— . 
The invariant as written in Eq. |2] is not sufficient to stabilize a first-order transition 
Pmmnl'— )'2221' which would have to occur at a single point in the parameter space. In 
order to stabilize this transition through a line in phase space, one needs to consider a 
strongly negative coupling term 7 and at least one positive six degree term (for example 
Spi). In such conditions, considering that the first coefficient ai follows the usual tempera- 
ture dependence ai=a(T — T^) with a > 0, and that rji is the driving order parameter, the 
critical behavior of 772 follows the triggering mechanism first discussed by Holakovsky^^ . 

The symmetry of this helical phase in such orthorhombic symmetry is extremely in- 
teresting. As already mentioned, the point group 2221', although not centrosymmetric, 
is not polar and therefore forbids a spontaneous polarization at the magnetic transition. 
Since the structure is incommensurate and the time-reversal operator is present in the point 
group, the linear magnetoelectric effect is also forbidden. However, invariants that are linear 
in the electric field (E) and quadratic in the magnetic field (H) are allowed by symmetry 
(whilst terms in HHE are forbiddden). The quadratic magnetoelectric effect is described 
by a third-rank tensor, fiijk symmetric in the last two indices and where i is the direction 
of the electric field and j and k indicate the direction of the magnetic fields. For the point 
group 2221' there are three terms allowed for this tensor that transform as the piezoelectric 
tensor, /3i23,/32i3 and /3312. One can easily show using the matrix representatives of A2 and 
A3 given in Tab. |T] that the following mixed invariants of the magnetic order parameters 
with E and H can be constructed: 

Fme = Pip2Sin{(t)) [PusE.HyH, + P2izEyH,H, + PsuE.H^Hy] (3) 
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FIG. 6. Phase diagram of the possible magnetic point groups that can be stabihzed from the para- 
magnetic group Pmmnl' by considering two order parameters ryi and 772 transforming respectively 
as the Ti and T2 representations (see text for details). 



This suggests that the application of low magnetic fields (if the system remains in the phase 
described by the symmetry of the zero-field phase) inclined with respect to the crystallo- 
graphic axes will induce an electrical polarization whose magnitude with be proportional to 
the square of the magnetic field. Of course, the application of sufficiently large magnetic 
fields could drive the system to other symmetry phases. Single crystal specimen are neces- 
sary to determine the behavior of the electric properties and magnetic properties under a 
magnetic field. 

It is important to realize that this result can be easily generalized to all incommensurate 
magnetic systems that possess a true magnetic handedness, i.e. all spirals in the unit-cell 
propagate with the same chitrality. In such systems and irrespective of the paramagnetic 
symmetry, all improper symmetry operations are lost in the magnetically ordered state and 
the time-reversal operator is preserved by the incommensurability, implying the at least one 
of the f3ijk is allowed and quadratic magnetoelectric effects are present. This magnetoelectric 
effect, even though of higher degree in the expansion of the free energy than the well-studied 
linear effect, is not bounded by thermodynamic conditions (see^ and references therein) 
and therefore interesting to investigate in prototype spiral magnets of the type discussed 
here. 
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In conclusion, measurements of the magnetic properties of Q;-CaCr204 and neutron 
powder diffraction have shown that this system is analog to the delafossite. Strong an- 
tiferromagnetic interactions in the triangular Cr02 planes leads to long-range magnetic 
ordering below Tjv=43K, with an incommensurate propagation vector k=(0,0. 3317(2), 0). 
The model for the magnetic structure, based on refinement of the neutron diffraction data 
and symmetry considerations, is consistent with an helicoidal modulation and a nearly 120° 
configuration in each triangular plaquette. The mirror plane symmetries are broken by the 
helicoidal arrangement, leaving the point group 2221', which allows for quadratic magneto- 
electric effects in all directions. 
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